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Abstract 

The evolution of surface morphology during the growth of N-polar (0001) GaN under N- 
rich condition is studied by kinetic Monte Carlo (kMC) simulations for two substrates mi scuts 
2° and 4°. The results are compared with experimentally observed surface morphologies of 
(OOOl) GaN layers grown by plasma-assisted molecular beam epitaxy. The proposed kMC 
two-component model of GaN(0001) surface where both types of atoms: nitrogen and gallium 
attach the surface and diffuse independently, explains that at relatively high rates of the step 
flow (miscut angle <2") the low diffusion of gallium adatoms causes surface instabilities and 
leads to experimentally observed roughening while for low rates of the step flow (miscut 4°), 
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smooth surface can be obtained. In the presence of almost immobile nitrogen atoms under 
N-rich conditions, the growth is realized by the process of two-dimensional island nucleation 
and coalescence. Additionally, we show that higher crystal miscut, lower crystal growth rate 
or higher temperature results in similar effect of the smoothening of the surface. We show that 
the surface also smoothens for the growth conditions with very high N-excess. The presence 
of large number of nitrogen atoms changes locally mobility of gallium atoms thus providing 
easier coalescence of separated island. 


Introduction 

The fact that in wurtzite N-polar GaN, the direction of the polarization fields is reversed compared 
to the Ga-polar direction, makes the N-polar GaN of interest for a variety of device applications 
such as high electron mobility transistors, solar cells, sensors and light emitting diodes.^i^ Identi¬ 
fication of the optimum growth conditions providing smooth surface morphology N-polar (0001) 
GaN is one of the key issues to enhance the development of these emerging applications. The 
growth of atomically flat N-polar (0001) GaN layers by plasma-assisted molecular beam epi¬ 
taxy (PAMBE) has been proven difficult due to the high adatom diffusion barriers because it is 
carried out at relatively low temperatures (750"C). Theoretical works showed that these diffu¬ 
sion barriers are substantially reduced when gallium adlayer is present on the surface during the 
growth.®^ However, previous experimental studies showed that the growth of the N-face GaN 
must be performed with less than one monolayer (ML) of excess Ga on the growing surface, in 
contrast to Ga-polarity, when more than 2ML can be stabilized. Therefore, a precise control of the 
growth conditions is needed to prevent the accumulation of Ga droplets.^ The morphologies re¬ 
ported for N-polar GaN layers grown by PAMBE under N-rich conditions on low miscut substrates 
were relatively rough either due to the hillocks or pattern formed by interlacing stripes covered by 
doubly-bunched atomic steps.^^In search for the optimum growth conditions of the device-quality 
N-polar GaN layers it is essential to understand of the role of surface miscut, growth temperature, 
growth rate and III/V ratio (N-excess) that may improve surface morphology. It is important to 
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note that in case of metal-organic vapor phase epitaxy (MOCVD), the N-polar (OOOl) GaN films 
for many years suffered from high surface roughness resulting from the formation of large hexag¬ 
onal hillocks The use of 4 deg miscut N-polar GaN substrates has proven to be effective method 
of improving surface morphology.!^ 

In this work the growth of N-polar GaN (0001) surface under N-rich conditions is investigated 
theoretically by kinetic Monte Carlo (kMC) simulations and compared to experimental results 
in order to identify the optimum growth conditions leading to smooth layer morphologies and 
understand the role of substrate miscut, temperature, growth rate and N-excess. We built two- 
component model of GaN (0001) surface. We use the model that is based on the one previously 
implemented to study growth and sublimation of GaN(OOOl) surface.^^^^^^The model is developed 
in such a way that both types of atoms building crystal: nitrogen and gallium are controlled. They 
attach to the surface, diffuse and detach independently. They can be adsorbed at the steps or form 
islands by nucleation at the terraces. In this work we study the influence of substrate miscut on 
the morphology of (0001) GaN layers by comparing the simulated systems to the morphologies of 
N-polar GaN 0001) layers grown bt PAMBE under N-rich conditions. 

In our kMC simulations of N-polar GaN (0001) growth under N-rich conditions we apply 
leV barrier for Ga surface diffusion as evaluated in ab-initio calculations® and 1.8eV high barrier 
for diffusion of N adatoms. With such a slow surface dynamics we model the crystal growth 
under large nitrogen flux at relatively high miscuts. We choose external fluxes and desorption 
rates in such a way that simulations follow experimental crystal growth rates. As a result the 
experimentally observed morphology patterns are well reproduced. We explain that the surface 
instability appears due to the too slow diffusion of gallium atoms at relatively high flux of incoming 
atoms and because of the wide terraces as compared to the diffusion length of Ga adatoms 
However if steps are located more densely, crystal surface becomes smoother under the same 
growth conditions. In two component system of GaN, at miscut of 4°, diffusion rate of N adatoms 
is below and of Ga adatoms close to the limiting value for the domain formation on the terrace.!^^^ 
If the diffusing particle has enough time to find a step edge before the next one attaches to it the 
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island nucleation process is suppressed. When this condition is fulfilled we observe smooth erystal 
growth realized either by step flow or mixed made of step flow and two dimensional nucleation 
process.^ It is interesting that we observe such growth mode for two component system, when 
one eomponent is slower and the seeond moves with rate elose to the expeeted limiting diffusion 
value. This limiting value ean be deereased by higher surfaee miseut, lower erystal growth rate or 
higher temperature. Our kMC simulations eonfirmed that all these methods lead to lower surfaee 
roughness. Thus we ean explain experimentally observed differences in surfaee morphologies of 
N-polar GaN (0001) grown at 2° and 4*^ surfaee miseuts. However, due to heterogeneous growth 
and immobility of nitrogen adatoms the eharaeter of stable and unstable surfaee evolution is a 
mixture of step flow and 2-D or 3-D island growth.^^^^® 

We find that low miseut angle (< 2 °) results in eharaeteristic surfaee instabilities (and sub¬ 
sequent high surface roughness), while for steeper cuts a smooth double-step surfaee pattern is 
observed. Other methods of surfaee smoothening lead to slightly different surfaee morphologies. 
For instanee higher temperature smoothens the surfaee, but eontrary to the high miseut ease no 
double step strueture ean be seen. Distanees between all steps are equal suggesting that at higher 
temperature all steps move with the same rate at low external partiele flux (slow erystal growth 
rate). The resulting surface morphology is very similar to this for 4° cut crystal. It can be seen 
when appropriate eorrelation functions are compared. However there are some differences in both 
patterns that suggest that dynamies of the proeess is slightly different. 

We have found that step roughness ean be smoothened also using increased nitrogen flux. As 
nitrogen praetically does not move along the surfaee at studied temperatures and the rate of gallium 
adatoms diffusion is also below limiting value, the character of this proeess is different than that 
studied above. The pattern that builds at the surfaee consists of islands that grow along step-edges 
and finally attaeh to the step. Resulting smooth erystal growth is realized by two dimensional 
island nueleation. We show that a flat surfaee is developed at relatively low miseut at high nitrogen 
rates. We diseuss numerieal data as a funetion of the growth temperature, growth rate defined by 
Ga flux and substrate miseut angle. The kMC results are compared with experimental data. 
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First in Section II the model is formulated, then in Section III we show how the surface structure 
changes with miscut, temperature, galium and nitrogen flux. Surface patterns obtained for the 
chosen parameters are compared to the experimental data. Finally on the basis of the surface 
dependence on the nitrogen flux we consider possibility of the surface smoothening via use of 
higher nitrogen fluxes. 

Kinetic Monte Carlo model 

Numerical simulations of the system were performed with use of two-component kMC model. It 
consists of gallium and nitrogen atoms arranged in a crystallographic lattice of N-face (0001) of 
gallium nitride. Two species interact via nearest neighbor (NN) and next nearest neighbor (NNN) 
forces which form N-N, Ga-Ga, and Ga-N bonds. Total energy of an adatom occupying one of the 
lattice sites is given by the sum over all NNs and NNNs surrounding it 

Ex=EGaN^ni + Exx ni ( 1 ) 

NN NNN 

where X corresponds to Ga or N atom and the value n/ depends on the state of the neighboring 
site. If the site is occupied ni = 1 otherwise n, = 0. Energies associated with interactions are 
Enn = 0-3, EcaGa = O.SSeV and EoaN = TbeV respectively. They were chosen in such a way that 
total energy of Ga atom built in the crystal is 10.6 eV. This value is in agreement with molecular 
dynamic calculations performed by Wang et. aliasing a bond order potential. 

Single simulation consists of several Monte Carlo (MC) steps. Each of them starts with par¬ 
ticle adsorption which occurs with probability equal to the external flux R Adsorbed particles are 
allowed to diffuse at the surface. That process is modeled in three stages. At first one of the neigh¬ 
boring sites is randomly chosen to be the target of a jumping particle. Next initial (E,) and final 
(Ef) energies of jumping atom are calculated using eq. (??). In the third stage of diffusion process 
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particles jump with probabilities: 


Pd = 


^d) when Ei>Ef 
ye^P^D otherwise 


where j8 = {kBT)~^ and v = sets the time seale. denotes the diffusion barrier that 

is different for X=Ga, N adatoms. It is known that diffusion barriers for both components under 
nitrogen rieh conditions on GaN(OOOl) surface are high.®^ We use barriers ealculated in Ref.,® 
and so Ga adatoms diffuse jumping over the barrier = leV and N adatoms jump over barrier 
B^ = l.SeV. Aeeording to the Eq. ?? partiele, whieh is strongly binded in its initial state, resides 
in a deep potential well and probability of the jump out is low. However when final state has low 
energy i.e. partiele tries to move to the deep well the energy barrier between these states deereases 
and jump rate is high. Finally jumps towards step from above oeeur with the same probability like 
those far from the step in sueh a way that effectively no additional barrier for jumps aeross step 
(Sehwoebel barrier) is present in the studied systems. During the last stage of a single MC step, 
the evaporation phenomenon is simulated. An assumption that ineoming atoms ean evaporate is 
neeessary in order to reaeh equal number of Ga and N atoms within the growing erystal. We set 
the following rate of particle desorption: 


Dx = ( 2 ) 

where dx is the desorption potential, different for Ga and N adatoms. We set dx — —0.7eV taking 
into aeeount that binding energy of single N atom at the surfaee Ei=5JeV while the same atom 
above the (0001) GaN surfaee eovered by nitrogen Ei = 0.9eV . This value of desorption barrier 
ensures that most of N adatoms that are not bonded to gallium atoms in the layer below, desorb 
from the surfaee, whereas the ones that have fallen onto the gallium layer stiek to the surfaee. For 
gallium we set dca = 0. Sueh ehoiee of desorption potentials provides that the number of adsorbed 
atoms of both types equalizes at the surfaee. After evaporation stage new MC step begins and the 
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above procedure starts from the beginning. 

At the beginning of each simulated process, the system consists of Ns straight steps. Helical 
boundary conditions in one direction and periodic in the second one are posed. Gallium nitride 
elementary cell consists of four mono-atomic layers, hence in order to close helical boundary 
condition properly. Ns is always divisible by 4. Additionally, because N-rich conditions are studied, 
we start the simulation with total N coverage of the surface, what changes during the growth 
process by evaporation of adatoms from the surface. 

Summarize, the kMC model used to describe N-polar GaN(OOOl) growth under N-rich condi¬ 
tions we build of two types of atoms: Ga adatoms that move slowly but are able to diffuse across 
terraces and reach steps, and N adatoms that are almost immobile, but can desorb from the sur¬ 
face and are frequently adsorbed because of the relatively high flux of N atoms. The consequence 
of such assumptions is that both atomic monolayers (nitrogen and gallium) can grow smoothly 
within some range of parameters, but due to totally different mechanisms. Gallium adatoms dif¬ 
fuse over the surface, and they start to build islands when incoming flux is too high, whereas 
nitrogen adatoms being practically immobile build islands all the time. Nitrogen islands have two 
different roles in the growth process: they build in the step thus completing the layer formation, 
but they also change diffusion of Ga adatoms. Below we analyze consequences of such adatom 
dynamics for the evolution of stepped surfaces. 


MBE experimental details 

In parallel we carried out the growth of 200nm thick N-polar GaN layers by PAMBE custom 
design Gen20A MBE reactor were active nitrogen was supplied from Veeco RE plasma source. We 
used bulk, commercially available GaN substrates from Saint Gobain with threading dislocation 
density around 5 x 10^ cm^^. Prior to growth, substrates were mechanically polished to get miscut 
angle of 2° and 4" toward [1010] direction and chemomechanically polished to obtain atomically 
smooth surface. Epi-ready substrates were then mounted by gallium to single 2 inch wafer close 


7 




Figure 1: (color online) (a) Surface pattern of N-polar GaN(OOOl) layer simulated by kMC on 2® 
miscut substrate, (b) Morphology of 200nm N-polar GaN layer grown by PAMBE under the condi¬ 
tion used in the kMC simulations on 2° miscut substrate. Fluxes of atoms are pQa = 4nmjmin and 
= I6nmlmin for both experiment and simulations. Temperature T = 150°C in both pictures. 



x[nm] 


Figure 2: (color online) (a) Surface pattern of N-polar GaN(OOOl) layer simulated by kMC on 4® 
miscut substrate, (b) Morphology of 200nm N-polar GaN layer grown by PAMBE under the condi¬ 
tion used in the kMC simulations on 4® miscut substrate. Fluxes of atoms are Foa = Anm/min and 
F^v = I6nm/min for both experiment and simulations. Temperature T = 150°C in both pictures. 
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to its center to ensure the same growth eonditions at both erystals. Growth of 200 nm thiek GaN 
was earried out at 750°C out under nitrogen rieh eonditions using nitrogen and gallium fluxes of 
Fn = I6nm/min and Foa = Anm/min (GaN equivalent growth rate) respeetively. After the growth 
samples morphology was examined by atomie foree mieroseope (AFM). 


Dependence of the surface ordering on the growth parameters 

Miscut 

[^presents the numerieal (a) and experimental (b) results of N-polar GaN layer morphology for the 
ease of 2° miseut surfaee. For both oases (experimental and oaloulations) growth parameters like 
temperature and atomie fluxes were the same. The size of ealeulated system is 105 nm x 105 nm, 
smaller than the area of AFM sean whieh is 1/im x 1/im. Despite this it ean be seen that in both 
eases orystal does not grow smoothly. In several eavities ean be easily seen, steps bend and 
beoome wavy leaving deep eavities on their way forward. In the prooess of further growth more 
profound struotures ean build around them. In[TJ) we ean see that the surfaee is rough, dominated 
by islands. This is typioal pattern expeoted for systems with slow adatom diffusion. 

The situation ehanges diametrieally when growth is performed under the same eonditions (im¬ 
pinging fluxes and temperature), but at the substrate eut 4°. presents the simulated surfaee 
pattern for 4° miseut. Double-steps with relatively straight edges are seen. presents the AFM 
sean of N-polar GaN layer grown by PAMBE. We ean see smooth, stepped surfaee. The double 
step strueture however is not so elearly visible like in the simulated plot above. The double step 
strueture is a eonsequenee of GaN lattiee geometry aeeording to whieh every seeond step perpen- 
dieular to [1010] direetion has different bond arrangement. Thus the rate of adatom adsorption at 
the step and resulting veloeity of every seeond step is different. 

In order to analyze and eompare surfaee parameters more preeisely we have ealeulated root 
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mean square (RMS) roughness 


RMS= (3) 

where h, is the z eoordinate (height) of i-th site over leveled surfaee and Hq is mean system height. 
RMS parameters for numerieal results presented in[^ and[^ are 0.38nm and 0.15nm respeetively. 
The same parameters ealeulated for eorresponding experimental data ([T]? and[^) are 6.11nm and 
0.49nm. Roughness measured for experimental surfaees is obviously higher due to the surfaee 
distortion or loeal defeets, however one ean notiee that RMS parameters for lower miseut are higher 
both in the experiment and simulations as well. Numerieally obtained surfaees are smoother also 
due to the shorter time of system evolution, whieh for experiments was about 40 minutes and 2-3 
minutes for simulations as well as due to different size of the surfaee area. 



Figure 3: (eolor online) Correlation funetions for experimental surfaees for a) 2° miseut and 
b) 4® miseut. The other parameters like in and Solid lines denote eorrelation ealeulated 
perpendieularly to to the initial step orientation - along x axis, and dashed lines denote eorrelations 
parallel to steps - along y axis. 


In order to illustrate and analyze surfaee features, in partieular their periodieity and amplitude 
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Figure 4: (color online) Correlation functions for simulated surfaces for a) 2° miscut and b) 
miscut. The other parameters like in[2and|^. Solid line denotes correlation calculated perpendic¬ 
ularly to to the initial step orientation - along x axis, and dashed line denotes correlations parallel 
to steps - along y axis. 


2° miscut 40 miscut 2° miscut + llx N flux 



x[nm] x[nm] x[nm] 


Figure 5: (color online) a) Simulated pattern of 2*^ miscut from Fig la. b) 4*^ miscut from Fig 2a 
and c) 2® miscut at high N flux Fn = 180nm/min and other parameters the same as in a). Red (dark 
gray) points denote nitrogen and green (light gray) gallium atoms. 
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we calculated correlation functions along and across steps according to the formula: 

Cin - fj) = YjiKri) - ho] [h{fj) - Hq] (4) 

i 

where the difference in position f/ — fj is calculated along x or y axis i.e. across (r = x) or along 
(r = y) steps, [^presents the correlation functions calculated for experimentally measured surface 
morphology of N-polar GaN(OOOl) layers grown by PAMBE on (a) 1° miscut and (b) 4® micut 
GaN substrate. We can see that the scale of variability of the 1° curve is much larger than that for 
4° surface i.e. ±30nm^ vs. ±0.2nm^ respectively. [^. shows that the correlation function calcu¬ 
lated for the surface of 2" oscillates in both directions. The length scale of its changes along and 
across steps is much longer than bilayer step distance (;^20nm). We can see that the characteristic 
correlation length is shorter across than along steps. The correlation function across steps of 4® 
miscut surface (|^) shows periodicity of around 10 nm length denoting bi-step surface structure. 
The function along step is flat and close to zero describing smooth surface. 

Correlation functions for simulated (0001) GaN surfaces on 2^^ miscut (j^) and 4° miscut (|^) 
exhibit many similarities to the experimentally corresponding counterparts. The characteristic 
length scales at 2° cut are shorter in both directions, thought it reflects less corrugated structure 
due to short available simulation times. Nevertheless the period of correlation function oscillations 
is still longer than bilayer step distance. Correlation function across the steps for miscut (0001), 
presented in|^, clearly shows the step structure. Oscillations of the function are very regular and 
all with the same amplitude what corresponds to the very regular structure. Thus the character of 
surface morphology is very well reflected by the shape of correlation function. 

The analysis of the correlation function for the (0001) GaN surface morphologies obtained 
on 2° miscut and A° miscut substrates confirms that the surface of a crystal grown at substrate 
with higher miscut is more regular then the surface which evolved at the substrate with lower 
miscut. In order to understand the causes of observed differences we carry out the following 
analysis: Formation of Ga 2D islands on the terraces is probable for 2° miscut because the average 
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diffusion length for Ga adatoms until they meet next partiele is around 8nm, whieh is elose to 
the terraee width for IP eut. To be more preeise we ealeulate the limit terraee length for the 
domain formation,that is given by the relation lcr= {D/py/^. For our data, due to the low 
diffusion barrier, the limit terraee length is around = 4a for nitrogen adatoms and ify — \1 a for 
gallium adatoms, where a denotes GaN lattiee eonstant a = 0.319nm, This gives = 1.3am and 

= 5.44am as eompared with the terraee width / = 8a at the surfaee of 4° miseut and Z = 16a 
in the 2® ease. This means that in this for 2° miseut (0001) GaN surfaee nitrogen adatoms have 
tendeney of ereating islands, whereas gallium adatoms are at the erossover length. However at the 
temperatures studied the surfaee has a tendeney to build bi-layer struetures. Henee terraee width 
that has to be eonsidered as being two times wider. As a result adatoms at miseut surfaee are 
already within the erossover region. 

In sueh situation, partieles adsorbed at the steeper surfaee start to form islands, but beeause 
these islands built up short distanee to step edges. Steps move forward due to the attaehment of 
other partieles. In sueh a way islands are absorbed at the steps . It ean be seen in |^. Green 
points that are plotted for gallium atoms are arranged in pattern that looks like many islands glued 
together. Nitrogen adatoms seem to have not mueh importanee here, beeause there are enough of 
them all around. Henee the system forms steps straight and distributed regularly in spaee. On the 
other hand, when miseut angle is low, terraees are wider and partieles need more time to reaeh the 
step. In sueh ease the probability of 2-dimensional nueleation is higher. As a result of this proeess 
islands would appear at terraees. They are partly adsorbed by steps, leaving some eavities beside 
them and irregular wavy pattern builds up ([^). Less mobile nitrogen adatoms that are adsorbed 
and desorbed seem to play a role in building 3D islands. On building next layer on top of island 
they ean be souree of effeetive Sehwoebel barrier for wandering gallium adatoms. For mueh higher 
nitrogen flux this reasoning does not work. Nitrogen adatoms are everywhere and they are rather 
smoothening the surfaee of growing 2° erystal as we ean see in|^. The surfaee morphologies 
of GaN (0001) layers grown on 2® miseut ([^ and on 4^^ miseut [^substrates elearly evolve in two 
direetions. From the mieroseopie point of view, there is a signifieant differenee between the time 
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Figure 6: (color online) (a) Surface pattern of N-polar GaN(OOOl) layer simulated by MC on 2® 
miscut substrate at high T. pQa = Anmlmm,FM = 16nm/min, T = SSO'^C . (b) Correlation functions 
calculated for this surface. Solid line denotes correlation calculated perpendicularly to to the initial 
step orientation - along x axis, and dashed line denotes correlations parallel to steps - along y axis. 
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Figure 7: (color online)) (a) Surface pattern of N-polar GaN(OOOl) layer simulated by kMC on 2° 
miscut substrate for low F. pQa = Inm/min, Fn = Anmimm, T = 150°C. (b) Correlation functions 
calculated for this surface. Solid line denotes correlation calculated perpendicularly to to the initial 
step orientation - along x axis, and dashed line denotes correlations parallel to steps - along y axis. 
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Figure 8: (color online) (a) Surface pattern of N-polar GaN(OOOl) layer simulated by kMC on 1° 
miscut substrate for high N flux. pQa = Anmimm, = 180nm/min, T = 150°C. (b) Correlation 
functions calculated for this surface. Solid line denotes correlation calculated perpendicularly to 
to the initial step orientation - along x axis, and dashed line denotes correlations parallel to steps - 
along y axis. 
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that adatoms need to reaeh atomie steps and the probability for eneountering one another. 

Temperature 

Higher miseut sueeessfully flattens the surfaee by shortening the time that adatom needs to reaeh 
the step. Knowing that we earned out another simulation for 2° miseut applying a higher growth 
temperature, thus inereasing surfaee diffusion rate to eheek if this will result in smoother surfaee. 
Inj^we show simulated surfaee morphology for (0001) GaN layer on 1° miseut surfaee obtained at 
the temperature 850°C, that is 100"C higher than the referenee (ef. Fig 1). In sueh ease diffusion is 
4 times faster than in the referenee example. It should eompensate 4 times longer lasting diffusive 
walk, what means effeetively 2 times longer terraee width. Comparing the surfaee morphologies 
in Fig. la and Fig. 6a we see the positive effeet of raised temperature on the surfaee smoothness. 
Caleulated RMS value for the N-polar (0001) GaN surfaee simulated at 850®C is 0.15nm, what is 
low value. However the morphologies obtained at raised temperature on 1° miseut (0001) substrate 
(|^ and at referenee temperature on 4° miseut (0001) GaN substrate are different. As ean be 
eoneluded from the eorrelation funetions, plotted in Fig. 4b and|^. Regular pattern aeross steps 
is elearly visible, but with no sign of double step strueture. In[^ the periodieity of the eorrelation 
funetion is around lOnm aeross the steps, while along the steps niee smooth terraees are indieated 
by low, elose to zero eorrelation line. Sueh strueture ean be also seen in Fig. 4b. However, when 
we look eloser at the surfaee morphology in Fig. |^we ean see that every seeond step has irregular 
strueture similar to that seen in Fig. [^. More detailed study shows that these island-like struetures 
are build by nitrogen and not by gallium atoms. Gallium steps at the other side have tendeney to 
bend at eharaeteristie length related to the diffusive distanee at given temperature. 

Growth rate 

The other method of surfaee smoothening is to use a lower flux F of ineoming partieles i.e. lower 
growth rate. In sueh a way diffusing adatom has more time to wander around on the surfaee and as 
a result the steps should beeome more straight. In|^ we see an effeet of the slow growth rate (flux 
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are 4 times slower than in the referenee surfaee in[^. Again we ean see a rather smooth surfaee. It 
differs from this obtained by inereased temperature. We ean see bi-step strueture, like in the ease 
of 4® miseut. Sueh strueture is refleeted in the shape of the eorrelation funetions presented in|^. 
Correlation in the direetion perpendieular to steps is very similar to this in|^ but in two times larger 
distanee seale. At the same time eorrelation along step quiekly deereases to zero. RMS = O.llnm 
is higher than in the ease of surfaee of erystal grown at higher temperature 


N-excess 

One more simulation was earried out to study the effeet of the inereased nitrogen flux on surfaee 
morpholog. In[^ we show the surfaee of GaN(OOOl) on 2° miseut obtained under the growth 
eonditions of high N exeess, i.e. the flux of ineoming nitrogen is 30 times higher than gallium 
flux. Part of this surfaee was already presented above in |^. It ean be seen that the surfaee is 
mueh smoother that this in however not in the same manner as A° miseut ease as seen in the 
simulated system in [^. It is easy to note that the inerease of nitrogen amount in the system leads 
to the improvement of surfaee quality. There are no eavities seen in and the height differenees 
are not so large as in the referenee (0001) GaN surfaee. However, the overall surfaee strueture is 
not so regular as in the eases studied before. The ealeulated value of RMS = 0.28nm is lower than 
the referenee 2° miseut GaN (0001) layer. When we look at the eorrelation funetion along steps 
presented in[^ we see that it deereases very quiekly, like for smooth struetures. The eorrelation 
funetion aeross the steps shows oseillations of eharaeteristie length for double step strueture i.e. 
20nm. In the last part of our study we analyzed the evolution of surfaee RMS parameter as a 
funetion of surfaee eoverage, eounted in the numbers of grown layers for three eases of (0001) 
GaN growth: high miseut, high temperature, and high N flux. The result is shown in Fig. For 
the ease of layer grown at the 4° miseut GaN substrate and high temperature 2° miseut growth 
proeess is evidently periodie. One ean see oseillations in RMS when the growth proeeeds. It ean 
be interpreted as layer by layer growth via ereation of islands, than absorption of these islands by 
approaehing steps, and the eyele ends up by next layer to be eompleted.^^^^^ We ean elearly see 
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Figure 9: (color online) RMS as a function of coverage shows periodicity of the process for the 
system of 4® cut and for high temperature growth for 2° cut compared to steady evolution for high 
nitrogen growth. 


periodicity of the length of one monolayer in the plotted curves. Typically, the higher temperature 
process is the smoother line we would expect, but it has the same period related to the number of 
layers. Third curve is plotted for high N-excess case RMS value stays at the same, relatively high 
level, showing that the growth process is realized in a different way. We attribute it to the island 
formation and coalescence that results in a stable RMS value over the growth time. 


Conclusions 

Kinetic Monte Carlo simulations have been employed to study the surface evolution of N-polar 
GaN(OOOl) layers grown on the substrates of 2® and A° miscut. The numerical results are com¬ 
pared with experimentally observed surface morphologies of GaN layers grown by MBE. Due to 
the high diffusion barriers during the growth of (0001) GaN under N-rich conditions the resulting 
surface is irregular and rough. Such surface can be observed for the crystal grown at 2° surface 
miscut at temperature around VSO'^C. Four different methods of smoothening the surface were in¬ 
vestigated using kinetic Monte Carlo method: (a) growth of GaN(OOOl) on high misut substrates, 
(b) at higher temperature, (c) at low growth rate and (d) using high N-excess. All of them lead to 
smoother crystal surfaces. The surface of A° cut seems to be the most regular one, that has been 
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also confirmed by by experimental results. Other approaches like: higher temperature or lower 
growth rate have also work as smoothening factor, but the character of the surface morphology 
is different than this for the high miscut surface. Their further evolution seems to happen less 
smoothly due to meandering steps and nucleation on terraces that lead to the formation of cavities. 

The surface that emerged during the growth at high N-excess using very high nitrogen flux is quite 
unusual, because microscopically it seems to be irregular, but correlation functions exhibit rather 
regular structure. Thus, we conclude that under high N-excess GaN(OOOl) growth happens regu¬ 
larly and the surface morphology remains smooth due to stable island formation and coalescence 
on the terraces. Therefore the use of high nitrogen flux for the growth of GaN(OOOl) under N-rich 
conditions seems to be a good method to obtain smooth surface morphology. 
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